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Abstract 

A closed expression for the density operator of the damped harmonic oscillator is 
extracted from the master equation based on the Lindblad theory for open quantum 
systems. The entropy and effective temperature of the system are subsequently calcu- 
lated and their temporal behaviour is surveyed by showing how these quantities relax 
to their equilibrium values. The entropy for a state characterized by a Wigner distri- 
bution function which is Gaussian in form is found to depend only on the variance of 
the distribution function. 



Shortened title: Damped quantum harmonic oscillator 
PACS numbers: 02.50.Ga, 05.20, 05.40.+j 



1 Introduction 



In the last two decades, the problem of dissipation in quantum mechanics, i.e. the consis- 
tent description of open quantum systems, was investigated by various authors [1-5] (for a 
recent review on open quantum systems see ref. It is commonly understood 0, 0] that 
dissipation in an open system results from microscopic reversible interactions between the 
observable system and the environment. Because dissipative processes imply irreversibility 
and, therefore, a preferred direction in time, it is generally thought that quantum dynamical 
semigroups are the basic tools to introduce dissipation in quantum mechanics. In the Markov 
approximation the most general form of the generators of such semigroups was given by Lind- 
blad This formalism has been studied for the case of damped harmonic oscillators ||^, 10 



and applied to various physical phenomena, for instance, the damping of collective modes in 



deep inelastic collisions in nuclear physics [11]. In [12| the Lindblad master equation for the 
harmonic oscillator was transformed into Fokker-Planck equations for quasiprobability dis- 
tributions and a comparative study was made for the Glauber P, antinormal ordering Q and 



Wigner W representations. In [13| the density matrix for the coherent state representation 
and the Wigner distribution function subject to different types of initial conditions were ob- 
tained for the damped harmonic oscillator. A remarkable feature resulting from the preceding 
analysis [n| is the agreement in form of the Lindblad master and Fokker-Planck equations 
with the corresponding equations familiar from quantum optics. Quantum optics is gener- 
ally formulated with the help of the quantum mechanics of the damped harmonic oscillator 
@, 1^ and the corresponding Brownian motion. Generally, the equation of motion treated 
in the theory of Brownian motion is the master equation satisfied by the density operator 
or the Fokker-Planck equation for the distribution function. The master and Fokker-Planck 
equations related to the problem of Brownian motion have been fully reviewed in the past 

i, %mM- 

In the present study we are also concerned with the observable system of a harmonic 
oscillator which interacts with an environment. The aim of this work is to explore further 
the physical aspects of the Fokker-Planck equation which is the c-number equivalent equa- 
tion to the master equation. Generally the master equation for the density operator gains 
considerably in clarity if it is represented in terms of the Wigner distribution function which 
satisfies the Fokker-Planck equation. We shall describe, within the Lindblad theory for open 
quantum systems, how the system under consideration evolves to a final state of equilibrium 
by calculating first an explicit form of the density operator satisfying the operator master 
equation based on the Lindblad dynamics. We subsequently derive the entropy and the ef- 
fective temperature of the quantum-mechanical system in a state characterized by a Wigner 
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distribution function which is Gaussian in form and show how these quantities relax to their 
equihbrium values. 

Entropy is a quantity which may be interpreted physically as a measure of the lack 
of knowledge of the system. The idea of effective temperature associated with the Bose 
occupation distribution has already been introduced |]T^ in relation to the entropy which 



has been derived in the framework of quantum theory of harmonic oscillator relaxation. The 
derivation of the entropy associated with an infinite coupled harmonic oscillator chain has 



similarly been elaborated for classical [|I^ and quantum mechanical systems |]I9| represented 
by a phase space distribution function. In the present paper we derive first the density 
operator of the damped harmonic oscillator in the Lindblad theory by applying a procedure 



similar to the techniques developed in the description of quantum relaxation [|T^, [T^, pT|] . 
When we denote by p{t) the density operator of the damped harmonic oscillator in the 
Schrodinger picture, the entropy S(t) is given by 

S{t) = —kTT{p In p) = ~k ^ < mlpln >< In p|m >, (1) 

m,n 

where k is Boltzmann's constant. While the explicit form of the general density matrix 



< m\p\n > is now available |]14[, the evaluation of the matrix element of the logarithmic 
operator In p is not an easy task. An alternative way of calculating the entropy is to com- 
pute straightway the expectation value of the logarithmic operator < In p >= Tr(plnp). 
Accordingly, the problem amounts to derive the explicit form of the density operator. 

The content of this paper is arranged as follows. In Sec. 2 we write the master equation 
for the density operator of the harmonic oscillator. Sec. 3 derives an explicit form of the 
density operator involved in the Lindblad master equation. Sec. 4 formulates the entropy 
and time-dependent temperature using the explicit form of the density operator and discusses 
their temporal behaviour. Finally, concluding remarks are given in Sec. 5. 

2 Master equation for the damped quantum harmonic 
oscillator 

The rigorous formulation for introducing the dissipation into a quantum mechanical system 
is that of quantum dynamical semigroups [|, |, 1). According to the axiomatic theory of 
Lindblad the usual von Neumann-Liouville equation ruling the time evolution of closed 
quantum systems is replaced in the case of open systems by the following equation for the 
density operator p: 

^ = i(t,(p)). (2) 
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Here, $t denotes the dynamical semigroup describing the irreversible time evolution of the 
open system in the Schrodinger representation and L the infinitesimal generator of the dy- 
namical semigroup $4. Using the structural theorem of Lindblad which gives the most 
general form of the bounded, completely dissipative Liouville operator L, we obtain the ex- 
plicit form of the most general time-homogeneous quantum mechanical Markovian master 
equation: 

where 

L{pit)) = -^[H,p{t)] + ^E([^.pw,^/] + [y^^pmi])- (4) 

Here H is the Hamiltonian of the system. The operators Vj and are bounded operators 
on the Hilbert space Ti. of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature 
are of this form after some rearrangement of terms, even for unbounded Liouville operators. 
In this connection we assume that the general form of the master equation given by (3), (4) 
is also valid for unbounded Liouville operators. 

In this paper we impose a simple condition to the operators H,Vj,Vj that they are 
functions of the basic observables q and p of the one-dimensional quantum mechanical system 
(with [q, p] = ihl, where I is the identity operator on Ti) of such kind that the obtained model 
is exactly solvable. A precise version for this last condition is that linear spaces spanned by 
first degree (respectively second degree) noncommutative polynomials in q and p are invariant 
to the action of the completely dissipative mapping L. This condition implies |^ that Vj are 
at most first degree polynomials in q and p and H is at most a second degree polynomial in 
q and p. Then the harmonic oscillator Hamiltonian H is chosen of the form 

H = Ho + ^(qp + pq), Ho = —p" + -^q\ (5) 
With these choices the Markovian master equation can be written [IIU 



dp XX z 

^ = -^[^o,p] - ^(A + /i)[g,pp + pp] + ^(A -/i)[p,pg + gp] 

-^[^, %P\] - [P, [P,P]] + [P,P]] + [P, %P]]l (6) 

where Dpp, Dqq and Dpq are the diffusion coefficients and A the friction constant. They satisfy 
the following fundamental constraints |T(]|: 

i)Dpp>0, ii)^,, >0, ill) DppDqq-Dlq>X^h'' /A. (7) 
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In the particular case when the asymptotic state is a Gibbs state 



Pg{oo) = e /Tre 



these coefficients reduce to 



^ X + Li, . huj ^ X — ah . fiuj ^ ^ 

Dpp = -^—Jimuj coth — — , = — coth — — , D„ = 0, 9 

2 2kl 2 muo 2ki 



where T is the temperature of the thermal bath. 

3 Evaluation of the density operator 

There are several techniques available to solve equations such as (6). By introducing the real 
variables xi^X2 corresponding to the operators 

jmuj 1 



m 



12, 13| we have transformed the operator form of the master equation into the following 



Fokker-Planck equation satisfied by the Wigner distribution function iy(xi,X2,t) : 

""^ E + l E Q'S^r^W, (11) 



i,j=l,2 ^-^^ ^ «J=1,2 'J^t'J^j 



where 



Since the drift coefficients are linear in the variables xi and X2 and the diffusion coefficients 
are constant with respect to xi and X2, Eq. (11) describes an Ornstein-Uhlenbeck process 



2^. Following the method developed by Wang and Uhlenbeck we solved |T3| this 
Fokker-Planck equation, subject to either the wave-packet type or the 5-function type of 
initial conditions. 

When the Fokker-Planck equation is subject to a Gaussian (wave-packet) type of the 
initial condition {xiq and X20 are the initial values of Xi and 0:2 at t = 0, respectively) 

Wy,{xi, X2, 0) = ^ exp{-2[(xi - xio)^ + (x2 - a;2o)^]}, (13) 
Tin 

the solution is found to be [ p^ : 

W^{xi, X2) = ■ ^^^rT- exp{-7^[0^(a;i - Xif + ^n>{x2 - ^2)^ + Xwi^l - Xi){x2 - X2)]},(14) 
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where 



Bn, = 9192 - ^93, 9i = 9*2 



^e^^* + ^(e^^*-l), ^73 = 2[e-^^* + ^(l-e-^^*)], (15) 

UJ A A 



9ia* +g2a - g^, ipn^ = 9i + 92 - 93, Xw = '2{9ia* + 92a) - 93{a + a*)- (16) 



We have put a = {fi — iQ)/uj, A = —X — iQ, di = {a^muDgg + 2aDpg + Dpp/muj) /Ti, 
d2 = {fTi^Dqq + 2fiDpq/uj + Dpp/muj) /% and Vt^ = uo'^ — /i^. The functions Xi and X2-, which 
are also oscillating functions, are given by [|l^ : 



Xi = e ^^[{xiQ (cos Vtt + ^ sin Vtt) + 0:20 77 sin Vtt] , 



(17) 



X2 = e ^^[{x2o{cosQt — ^sinQt) — xiQ—sinQt]. (18) 

The solution (14) of the Fokker-Planck equation (11), subject to the wave-packet type of 
initial condition (13) can be written in terms of the coordinate and momentum (< A >= 
Tr(pA) denotes the expectation value of an operator A) as |1T0|] : 



1 1 

W{q,p) = exp{- — [0(g- <q>Y + tpip- <V>? - 2x(g- < q >)(p- < V >)]},(19) 



where 



<q>= e"^* [(cos fit + ^sin^]^) < g(0) > H — ^ sin fit < ]3(0) > 



n 



(20) 



<P>=e —SI 



sin fit < g(0) > + (cos ~ ^ sin fit) < p(0) >] , 



(21) 



= cr, 



pp 



--< > - < p >^= 



■^fi^*^"" 



(22) 



ip = 



qq 



--< <f > - < q >^= 



(23) 



1 



X = (^pqit) = -< qp+pq > - < q>< p>= -^Xw, 



-X 



(24) 



To get the explicit expression of the density operator, we use the relation p = 27ihN{Ws{q, p)}, 
where Wg is the Wigner distribution function in the form of standard rule of association and 
A'" is the normal ordering operator [^, ^ which acting on the function Ws{q,p) moves 



5 



all p to the right of the q. By the standard rule of association is meant the correspondence 
pm^n _^ ^npm i^g^^ggj^ functions of two classical variables {q,p) and functions of two quantum 
mechanical canonical operators {q,p). The calculation of the density operator is then reduced 
to a problem of transformation of the Wigner distribution function by the N operator, pro- 
vided that Ws is known. A special care is necessary for the N operation when the Wigner 
function is in the exponential form of a second order polynomial of q and p. The Wigner 
distribution function (19) corresponds however to the form of the Weyl rule of association 



25 1 . This function can be transformed into the form of standard rule of association by 



Ws{q,p) = exp{Uh^)Wiq,p). (25) 

Upon performing the operation on the right-hand side, we get the Wigner distribution func- 
tion Ws, which has the same form as the original W multiplied by h but with x — 'i/;,/2 in 
place of X '■ 

^'^^'^^ " 2^^^^^"^^*^^^" < ^ >)^ + ^(P- <P>?- 2x'(g- < q >){P- < P >)]},(26) 
where 

^ = 0^ - x'^ X =X~ i^- (27) 
The normal ordering operation of the Wigner function Wg in Gaussian form can be carried 



out by applying McCoy's theorem |2^, |2^ , which states that 



[Jexp(-z/;,7)]^/^exp(ag^ + (3p^ + -fqp) = A^[exp(Ag^ + 5/ + Gqp)], (28) 

where a = A/C, f3 = B/C, C = sinhT/rj, T = -ih{-f^ - 4a(3y/\ with J = coshT + 
i^7sinhr/r = 1/(1 — ihG). Having performed a straightforward calculation, we get the 
explicit form of the density operator: 



n ,1, ^ 1 . , ^ 



2 



^2 ^-zhx' 2hJi-inx' + W ^i^-^^x') 



h 

x[0(g- < q >f + %lj{p- < p >f - (x' + «-)[2(g- < q >){p- <p>)-ih]]}, (29) 

The density operator (29) is in a Gaussian form, as was expected from the initial form of the 
Wigner distribution function. While the density operator is expressed in terms of operators 
q and p, the Wigner distribution is a function of real variables q and p. When time t goes to 
infinity, the density operator approaches to 
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P(oo) = / ^ , exp{-— ^ln ^^^|| [ffpp(oo)g^ + agg(oo)p^-ffpg(oo)(gp + j3g)]}, (30) 
where a = (Tpp(oo)(Tgg(oo) — dp (oo) and [llO| : 



1 



oo 



"^■^'^^^ 2(mcu)2A(A2 + r]2) 



[(mcu)2(2A(A + /i) + cu2)D,, + iu'^Dpp + 2mcj2(A + fx)Dpg], (31) 



app(cx)) = 2A(A2V Q"^) [("^^^^^^^gg + (^'^('^ - Ai) + ^^)Dpp - 2mu;^(A - /i)L'pg], (32) 



(^pgioo) = ^-^^^^^—^[-(\ + fx){mufD,, + (A - /i)Z)pp + 2m(A2 - fi^)Dpg]. (33) 

In the particular case (9), 

, , h ^ hu , , hmuj Tiuj / \ r. /^A\ 
foo(oo) = coth— — , ariJoo) = coth— — -, crr,o(oo) = (34) 

and the asymptotic state is a Gibbs state (8): 

Pg(oo) = 2smh^exp{- — (— p + )}. (35) 

4 Entropy and effective temperature 

Because of the presence of the exponential form in the density operator, the construction of 
the logarithmic density is straightforward. In view of the relations (22)-(24), the expectation 
value of the logarithmic density becomes 

< Inp >= In/i - -In 5 ^ In — -. 36 

By putting hv = \/5 — h/2, we finally get the entropy in a closed form, 

S{t) = k[{u+l)\n{u + l)-u\nu]. (37) 

Because of the identity S = function z/ takes explicitly the form 



where 



= exp(-4At)(2^Re4 ^ + 111 _ ^ + 2^, 



to A |A|2 A2 A 

-2exp(-2At)[Re ( + ^) exp2^^]^ ) - ^ + ^] + ^ - I- (39) 



It is worth noting that the entropy depends only upon the variance of the Wigner distribution. 
When time t 00, the function u goes to s = uj{dyX^ - \di\'^/{\^ + fi2))i/2/2Q _ 1/2 and 
the entropy relaxes to its equilibrium value S{oo) = k[{s + 1) ln(s + 1) — sins]. 

A similar expression to (37) has been formerly obtained |]T^ |18| in the framework of 
the quantum theory of oscillator relaxation. It should also be noted that the expression (37) 
has the same form as the entropy of a system of harmonic oscillators in thermal equilibrium. 
In the later case u represents, of course, the average of the number operator [|TU|. Eq. (37) 
together with the function u defined by (38), (39) is the desired entropy for the system. While 
the formal expression (37) for the entropy has a well-known appearance, the explicit form of 
the function u displays clearly a specific feature of the present entropy. We see that the time 
dependence of the entropy is represented by the damping factors exp(— 4At), exp(— 2At) and 
also by the oscillating function exp 2iQt. The complex oscillating factor exp 2iQt occurs also 



in the second moments, as shown in [|l^. This factor reduces, however, to a function of the 
frequency u, expressly exp2iut for /x — > or if ^ 1 (i.e. the frequency uj is very large 
as compared to /i). 

In the case of a thermal bath (8), (9), we may define a time-dependent effective temper- 
ature Te, by remarking that at infinity of time the quantity u goes, according to (2.22) in 
ref. [0, to the average thermal phonon number < n >= {exp{hu! / kT) — 1)^^. Thus u may 
be considered as the time variation of the thermal phonon number. Accordingly we may put 
in this case 

(exp^-l)-^ = ^. (40) 

The function u, which goes to < n > as t tends to infinity, vanishes at t = 0. From (40) the 
effective temperature can be extracted as 

/c[ln(z/ + 1) — In z/J 

In terms of effective temperature we may say that the system at time t is found in thermal 
equilibrium at temperature Tg. Then, in terms of the temperature, the entropy takes the form 
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This form of the entropy for the thermally excited oscillator state is well-known. The effec- 
tive temperature approaches thermal equilibrium with the bath, Tg ^ T, as the value of t 
increases. 

5 Concluding remarks 

Recently we assist to a revival of interest in quantum brownian motion as a paradigm of 
quantum open systems. There are many motivations. The possibility of preparing systems 
in macroscopic quantum states led to the problems of dissipation in tunneling and of loss 
of quantum coherence (decoherence) . These problems are intimately related to the issue of 
quantum-to-classical transition. All of them point the necessity of a better understanding 
of open quantum systems and all requires the extension of the model of quantum brownian 
motion. Our results allow such extensions. The Lindblad theory provides a self consistent 
treatment of damping as a possible extension of quantum mechanics to open systems. In 
the present paper we have studied the one-dimensional harmonic oscillator with dissipation 
within the framework of this theory. We have first obtained the density operator from the 
master and Fokker-Planck equations. The density operator in a Gaussian form is a function 
of the position and momentum operators in addition to several time dependent factors. The 
explicit form of the density operator has been subsequently used to calculate the entropy and 
the effective temperature. The temporal behaviour of these quantities displays how they relax 
to the equilibrium value. In a future work we plan to use the entropy of the system in order 
to quantify the extent of decoherence in the course of the system-environment interaction. 
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